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Abstract. In this paper, we consider the three-dimensional inhomogeneous 
Navier-Stokes equations with density-dependent viscosity in presence of vac¬ 
uum over bounded domains. Global-in-time unique strong solution is proved 
to exist when ||Vmo||l 2 is suitably small with arbitrary large initial density. 
This generalizes all the previous results even for the constant viscosity. 
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1. Introduction 


The Navier-Stokes equations are usually used to describe the motion of fluids. 
In particular, for the study of multiphase fluids without surface tension, the 
following density-dependent Navier-Stokes equations acts as a model on some 
bounded domain 12 C R N (N = 2, 3), 


( 1 . 1 ) 


' p t + div (pu) = 0, in 12 x (0, T ], 

(pu)t + div (pu ® u) — div (2 p(p)d) + VP = 0, in 12 x (0, T], 
< div u = 0, in 12x[0,T], 
u — 0, on (912 x [0, T], 


k p\t=o — Po, u\t=o — Wo, in 12. 

Here p, u, and P denote the density, velocity and pressure of the fluid, respectively. 


d = i [Wu + (Vuf] 


is the deformation tensor. 


p = p(p) states the viscosity and is a function of p, which is assumed to satisfy 
(1.2) p e C^O, oo), and p > p > 0 on [0, oo) for some positive constant p. 
In this paper, we study the initial boundary value problem to the system (II.lj) - 

(L2D ■ 

The mathematical study for nonhomogeneous incompressible flow was initiated 
by the Russian school. They studied the case that p(p) is a constant and the 
initial density p 0 is bounded away from 0. In the absence of vacuum, global 
existence of weak solutions as well as local strong solution was established by 
Kazhikov mm- The uniqueness of local strong solutions was first established 
by Ladyzhenskaya-Solonnikov [22] for the initial boundary value problem, see 
also [25]. Furthermore, unique local strong solution is proved to be global in 

l 








2 


XIANGDI HUANG AND YUN WANG 


2D [26]. In recent years, Danchin initiated the studies for solutions in critical 
spaces. He [51 ITU] derived the global well-posedness for small initial velocity in 
critical spaces, where density is close to a constant. For some subsequent works, 
refer to PM and references therein. We remark that in the very interesting 
papers mum Danchin-Mucha studied the case for which density is piecewise 
constant, see also some generalizations in 2D [IT8] . 

When initial vacuum is taken into account and /i(p) is still a constant, Simon 
[26] proved the global existence of weak solutions. Later, Choe-Kim [7] proposed 
a compatibility condition as (II.4p below to establish local existence of strong 
solution. Global strong solution allowing vacuum in 2D was recently derived 
by the authors [20]. Meanwhile, some global solutions in 3D with small critical 
norms have been constructed, refer to the results in [2[|8] and references therein. 

Finally, we come to the most general case: viscosity p(p) depends on density p. 
Most results were concentrated on 2D case. Global weak solutions were derived 
by the revolutionary work [Hi, 23] of DiPerna and Lions. Later, Desjardins [13] 
proved the global weak solution with more regularity for the two-dimensional 
case provided that the viscosity function p(p) is a small pertubation of a positive 
constant in L°°- norm. Very recently, Abidi-Zhang [3] generalized this 2D result 
to strong solutions. Regarding the 3D case, Cho-Kim [6] constructed a unique 
local strong solution by imposing some initial compatibility condition. Their 
result is stated as follows: 

Theorem 1.1. Assume that the initial data ( po,uo ) satisfies the regularity con¬ 
dition 

(1.3) 0 < p 0 G W 1 ’ 9 , 3 < q < oo, u 0 G Hl a D H 2 , 
and the compatibility condition 

(1.4) - div (p,(p 0 ) [Vu 0 + (Vu 0 ) T ]) + VP 0 = p|p , 

for some (Pq, g ) G H 1 x L?. Then there exists a small time T and a unique strong 
solution (p, u, P) to the initial boundary value problem (II. ip such that 

peC([0,T]-, IF 1 ’ 9 ), V«, P G C([0, T]; if 1 ) n L 2 (0,T; W 1,r ), 

Pt e C([0, T]; L q ), y/put G L°°(0, T; L 2 ), u t G L 2 (0, T; if 1 ), 

for any r with 1 < r < q. Furthermore, if T* is the maximal existence time of 
the local strong solution (p, u ) , then either T* = oo or 

(1.5) sup (||Vp(f)||L? + ||Vu(f)|| L 2 ) = oo. 

0<t<T* 

Motivated by the global existence result [8j for the special case that p is a con¬ 
stant, we aim to establish global well-posedness result for variable coefficient case. 
However, due to the strong coupling between viscosity coefficient and density, it’s 
more complicated and involved with variable coefficient p(p) and requires more 
delicate analysis. 

Our main result proves the existence of global strong solution, provided || X7u 0 \\ L 2 
is suitably small. To conclude, we arrive at 
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Theorem 1.2. Assume that the initial data ( Po,u 0 ) satisfies (II.31) - (II.4p . and 
0 < Po < P- Then there exists some small positive constant eo, depending on 12, 
q,p,~p = sup n(p), p, || V/t(p 0 ) ||x,® (= M), such that if 
[ 0 , 75 ] 

(1.6) ||Vii 0 ||l 2 < eo, 

then the initial boundary value problem ijnj admits a unique global strong solu¬ 
tion ( p,u ), with 

p G C([0, oo); W 1 ’*), Vm, P G <7([0, oo); H 1 ) n L 2 oc ( 0, oo; H^), 
p t G C([0, oo); L q ), J-pu t G L% c ( 0, oo; L 2 ), u t G L 2 loc ( 0, oo; Hi), 

for any r, 1 < r < q. 

The main idea is to combine techniques developed by the authors in [8lfT9] and 
time weighted energy estimates successfully applied to compressible Navier-Stokes 
equation by Hoff im 

Let’s briefly sketch the proof. First we assume that ||V//(p)||l<! is less than 4 M 
and ||Vu|| 2 2 is less than 4||Vu 0 || 2 2 on [0, T\, then we prove that in fact ||V/x(p)||l 9 
is less than 2 M and || Vu|| 2 2 is less than 2|| Vuoll^ on [0, T], under the assumption 
||Vu 0 ||l 2 < e o < \- On the other hand, the control of ||V//(p)||l« and ||V«||l 2 
lead to uniform estimates for other quantities, which guarantee the extension of 
local strong solutions. All the above procedures depends on a time independent 
bound of ||VuH^i^oo. Thanks to bounded domain, u indeed has exponentially 
decay rather than insufficient polynomial decay for the whole space, that’s the 
main reason why we can only treat system (II. ip in bounded domain. 

Remark 1.1. Theorem ! 1 . 2\ also hold for 2D case. Since the proof is quite similar, 
we omit it for simplicity. 

The rest of the paper is organized as follows: Section 2 consists of some nota¬ 
tions, definitions, and basic lemmas. Section 3 is devoted to the proof of Theorem 

o 


2. Preliminaries 

12 is a bounded smooth domain in M 3 . Denote 

[ f dx = f f dx. 

J Jn 

For 1 < r < oo and k G N, the Sobolev spaces are defined in a standard way, 

L r = L'{Sl), W k ' = {/ e V : V‘/£l r ), 

H“ = W kJ . = {/ 6 C 0 “ : div / = 0}. 

Hi = Cq°, Hl a = Cjff, closure in the norm of H 1 . 

High-order a priori estimates rely on the following regularity results for density- 
dependent Stokes equations. 
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Lemma 2.1. Assume that p G W 1,q , 3 < q < oo, and 0 < p < p. Let ( u,P ) G 
H L x L 2 be the unique weak solution to the boundary value problem 



(2.8) — div (2 p(p)d) + VP = F, div u — 0 in 12, and 


where d = \ [V« + (Vu) T ] and 


/x G C 1 [0, oo), p < p(p) < p on{0,p\. 


Then we have the following regularity results: 
(1) If Fe L 2 , then ( u,P ) G H 2 x H 1 and 



(2.9) 


where d 2 satisfies 




(2) If F G L r for some r G (2 ,q), then ( u,P ) G W 2,r x W 1,r and 



( 2 . 10 ) 


where 



r 


1 

<? 


Pere t/ie constant C in (I2.9j) and fl2.10p depends on 12, q, r. 

The proof of Lemma 12.11 has been given in [6j , although the lemma is slightly 
different from the version in [6]. We sketch it here for completeness. 

Proof. For the existence and uniqueness of the solution, please refer to Giaquinta- 
Modica [16]. We give the a priori estimates here. Assume that F G L 2 . Multiply 
the first equation of (12. 8 p by u and integrate over 12, then by Poincare’s inequality, 


J 2p(p)\d\ 2 dx 


I 


F ■ udx < ||P|| i2 • \\u \\ L 2 < C\\F \\ L 2 ■ ||Vu|| i2 . 


Note that 



hence 


( 2 . 11 ) 


IIVuiu* < CV'IIOA 


Since 




0, according to Bovosgii’s theory, there exists a function 


v G Hq, such that 
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and 


||Vn|| L2 < C 


P 


p(p) 


L 2 


Multiply the first equation of (j2,8[) by —v, and integrate over 0, then 
f P 2 r r 


rip) 


dx = — J F ■ v dx + 2 J n(p)d '■ V v dx 
< \\F\\l 2 ■ Ml 2 + Cjl • || V?i|| L 2 • || Vn|| L 2 

<C'||F|| L 2-||Vn|| L 2 + C'^-||F|| L 2-||Vn|| L 2 

P 


< C£ ■ IIFIItf • 
P 


On the other hand side, 


Hence, 

( 2 . 12 ) 


P 


P{p) 


L 2 


P 2 , f P 2 , 
dx > n I , dx. 


P{p) 


P 


p(py 


p(p) 


<c4 

L 2 P 


Ill’ll 


L 2 ■ 


The first equation of (12.8ft can be re-written as 


— A u + V 


P 


d 


rip)J p(p) 


F 2d-V^(p) PV rip) 


p{p) 


p{py 


By virture of the classical theory for Stokes equations and Gagliardo-Nirenberg 
inequality, we have 


\\ u \\ H 2 + 

< C 


V 

F 


(- 


P 


P{p) 


\P(P) 

+ 


L 2 


L 2 

d ■ Vrip) 


p{p) 


+ 


L 2 


PV^(p) 


a(p) : 


L 2 

-1 


< C p^\\F\\ L 2 + fi-'WVnl \l, ■ ||Vm|| J, + p-'WVp || i9 

V- - ~ 

< C [ff'WF || L 2 +^- 1 ||V A i||m • ||Vn||^ 2 2 • ||V«||^ 2 

+ p Iv/iIIm • 


2q 

L 9-2 


p 

0 2 

„ fP\ 

1 @2 

— 


V — 


p 

L 2 

V p J 

L 2 
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By Young’s inequality, 


Ml #2 + 


(2.13) 


—) 

mp)J 


L 2 


<C^\\F\\ L , + C^\\S7^\\ e L \-{ ||Vu||l» + 


L 2 


< CyT'WF |M + CpT^ 1 • ( 1 + h - ) • IIVM1 2 , • \\F\\l 2 


< C I \ + ^TillV^llS ] ||F||y, 


where 6b satisfies 


Similarly, 


1 1 _ 1 n _ Q 

2 q 3 6’ q- 3 


(2.i4) n«in a , + ii v (p/M)\\ Lr < c | i + -^hvmiis i iifiu-, 


/i jjB r 


where 


5 _ 1 

6 r 

1 _ 1 
3 q 


□ 


3. Proof of Theorem 11.21 

The proof of Theorem 11.21 is composed of two parts. The first part contains 
a priori time-weighted estimates of different levels. Upon these estimates, the 
second part uses a contradiction induction process to extend the local strong 
solution. The two parts are presented in Subsections 3.1 and 3.2, respectively. 

3.1. A Priori Estimates. In this subsection, we establish some a priori time- 
weighted estimates. The initial velocity belongs to H 1 , but some uniform esti¬ 
mates of higher order and independent of time are required. To achieve that, we 
take some power of time as a weight. The idea is based on the parabolic prop¬ 
erty of the system. In this subsection, the constant C will denote some positive 
constant which maybe dependent on U, q, but is independent of p 0 or uq. 

First, as the density satisfies the transport equation (11.lib and making use of 
ms, one has the following lemma. 

Lemma 3.1. Suppose (p,u,P) is the unique local strong solution to (11.ip on 
[0, T\, with the initial data (p 0 ,u 0 ), it holds that 

0 < p(x,t) < p, for every (x,t) G U x [0, T], 

Next, the basic energy inequality of the system (11,11) reads 
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Theorem 3.2. Suppose ( p,u,P ) is the unique local strong solution to (11.11) on 
[0,T], with the initial data (p 0 ,u 0 ), it holds that 

(3.15) f p\u(t)\ 2 dx+ f f p(p)\d\ 2 dxds < C-~p- ||m 0 ||l 2 , for every ie[0,T], 


or in other words, 


(3.16) / p\u[t)\ 2 dx + p 



|Vw| 2 dxds < C■ p■ \\u 0 \\ 2 L 2 , for every t G [0, T]. 


Proof. The proof is standard. Multiplying the momentum equation by u and 
integrating over D yield that 

1 d 


2 dt 


p\ufdx + 2 / p{p)\d\ 2 dx = 0 . 


Then (13.161) is true owing to the fact 2 J \d\ 2 dx = J \Wu\ 2 dx and p(p ) > p. □ 


Denote 


and 


M 2 

M r 


M = ||V/i(p 0 )|| 

1 p 1 

P + p p l / d 2+i 

1 p 1 

p + p p l / dr+l 


Li, 

■ (4 M)^, 



Theorem 3.3. Suppose ( p,u,P ) is the unique local strong solution to (II. ip on 
[0,T], with the initial data (po,uo), and satisfies 


(3.17) 

sup \\X7p(p(t))\\ Lq < 4M, 
te[o,T] 

and 


(3.18) 

sup ||Vu(t)||ia < 4 ||Vh 0 || 2 L 2 < 1. 
te[o,T] 

There exists 

a positive number C\, depending on D, q such that if 

(3.19) 

C l p~ 2 (Ml + ■ \\\7u 0 \\ 2 L 2 < In 2, 

then 


(3.20) 

- [ WVpUtWh dt + sup ||Vu(t)||i a <2||V« 0 ||i a . 

T Jo te[o,T] 


Before the proof of Theorem 13.31 let us introduce an auxiliary lemma, which is 
a result of the iy 2,2 -estimates in Lemma 12.11 

Lemma 3.4. Suppose (p, u, P ) is the unique local strong solution to (11.11) on 
[0, T] and satisfies 

sup ||V/i(p(t))|| L9 < AM. 
te[o ,T] 

Then it holds that 


||Vm|| h i < CM 2 \\p Ut \\ L i + CM 2 ■ p 2 ■ || Vn||l 2 . 
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Proof. The momentum equation can be rewritten as follows, 

(3.21) — 2div (p(p)d) + VP = — pu t — (pu ■ V)w. 

It follows from Lemma [2.11 and Gagliardo-Nirenberg inequality that 

||Vtt||jfi < CM 2 (llpMilU 2 + II pu ■ Vu|| lO 

< CM 2 \\pu t \\ L 2 + CM 2 ■ p • ||m||l 6 • || Vw|| L 3 

< CM 2 \\pu t \\ L 2 + CM 2 -p ■ \\Vu\\l 2 • UVulll,. 
By Young’s inequality, 

IIVm||hi < CM 2 \\pu t \\ L 2 + CM 2 ■ ~p 2 • ||Vu|| 3 i2 . 


Proof of Theorem \‘3.S\ Multiply the momentum equation by u t 
f1, then 

[ p\u t \ 2 dx + 






< 


dt „ 

pu ■ X7u ■ u t dx 


p(p)\d\ 2 dx 

+ C [ |V/i(p)| • |«| • | Vu| 2 dx. 


Here we used the renormalized mass equation for p(p), 

d t [p(p)} + u-X7p(p) = 0, 

which is derived due to the fact cliv u = 0. 

Applying Gagliardo-Nirenberg inequality and Lemma [3.41 

^ pu ■ X7u ■ u t dx 

< \\\-MtWh+CT>- IIVulll,' ||Vu||ai 

< gllvMh + CJ. HVulli. ■ [CM 2 \\fm t \\ L , + CM 2 V||V«||| 2 ] , 

and similarly, 

C J |V/r(p)| ■ \u\ ■ | Vu| 2 dx 
<C\\Vp(p)\\ L3 -\\u\\ L6 -\\Vu\\ 2 L4 

< C\\Vfi(p)\\ L 3 ■ ||Vu||f a ' IIVulll, 

< CM\\Vu\\l 2 • [CM 2 \\pu t \\ L 2 + CM 2 p 2 \\Wu\\ 3 l2 \K 
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Hence, by Young’s inequality, 

d 


J p\ u t\ dx + — J p{p)\d\ 2 dx 
< IWsfrhfv + \\\-Mt\\h + CMlf\\xjuf L2 + \\W-mWh 


C (MMj ■ pi||Vu||| a Y + CMMl ■ f ■ ||V«||® a 


+ 


< x \WpMI> + C {Ml + M 4 M 2 6 + MM 3 ) • p 3 • ||V U ||« a . 


So we have 


(3.22) 


d 


p\u t \ 2 dx + - / p{p)\d\ 2 dx < C {Ml + M 4 M 2 6 ) • p 3 • \\Vu\\ 6 l2 . 


Integrate with respect to time on [0, t], 
1 rt 
P Jo 

Applying Gronwall’s inequality, 

1 



p\u t \ 2 dxds+ sup / |Vn| 2 dx < Cp' 1 {Ml + M 4 M 2 6 ) -p 3 • / || Vu||| 2 ds. 
se[o,t] J Jo 


P Jo 


p\u t \ dxdt + sup / \S7u\ 2 dx 
te[o,T] J 


< || VmqII^ • exp l Cp~ 1 {M'l + M 4 M 2 6 ) • p 3 


II Vm||^ 2 dt 


According to Theorem 13.21 and the assumption (13.18(1 . 


(3.23) 


Hence, we arrive at 


(3.24) 


P Jo 


IIVu||^ 2 dt < sup ||Vm || 2 2 • / ||Vm|| 2 2 ds 

te[o,T] Jo 

< Cp - 1 - P- |K ||£ 2 

< C/Y 1 - p- 11V ti 0 11 z, 2 - 


p\u t \ 2 dxdt + sup / \S7u\ 2 dx 

tEfOjT 1 ] J 


< II'Vtx 0 ||| 2 exp{C 1 p- 2 (M 2 2 + M 4 Mf) • p 4 • ||V« 0 ||i 2 }. 

Now it is clear that (I3.20p holds, provided (I3.19j) holds. 

As a byproduct of the estimates in the proof, we have the following result. 


□ 


Theorem 3.5. Suppose ( p,u,P ) is the unique local strong solution to (11.11) on 
[0,T], with the initial data {po,uo), and satisfies the assumptions (13.171) - (13 .1 9h as 
in Theorem \3.3l Then 

rT 


(3.25) 


P Jo 


C ■ p 


t\\y/pU t \\ L 2 dt + sup t||V«||i 2 < 

te[o,T] p 


|| Vu 


o|Il2. 
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Proof. Multiplying (j3.22[) by t, as shown in the last proof, one has 


(3.26) 


I 1 Jo 


< 


P Jo 

According to Theorem 13.21 


t\\y/pu t \\ 2 L 2 dt + snp £||Vm|| 2 2 
te[o,T] 

rp 

p(p)\d\ 2 dxdt ■ exp{Ci/i~ 2 (Mf + M 4 M 2 ) • p 4 • || Vw 0 |Il 2 }- 


(3.27) 

Hence, 


p(p)\d\ 2 dxdt <C ■ p- \\uq\\ 2 L 2 < C ■ p ■ || Vn 0 || 2 2 . 


P Jo 


tWy/putWtf dt + sup t||Vn||| 2 
te[o,r] 


< — ||Vu 0 ||i 2 ■ exp{C 1 /i- 2 (M| + M 2 M 2 6 ) • p 4 • ||Vw 0 ||| 2 } 

P 

< ^||Vix, 


P 


2 

0||Z,2. 


□ 


Theorem 3.6. Suppose ( p,u,P ) is the unique local strong solution to (11.11) on 
[0,T], with the initial data (p 0 ,u 0 ), and satisfies the assumptions (13. 17|) - 03.191) . 
Then 

r r-T 


(3.28) 

and 


sup t / p\u t \ 2 dx + p / t||Vn t || 2 2 


te[o,T] 


<C'||Vno||l 2 e 1 -exp{C'e 2 } 


(3.29) 

where 

(3.30) 


©i = 


sup t 2 / p\u t \ 2 dx + p / t 2 || 'Vu t \\ 2 L 2 

t 2 G[0,T] J ~J 0 

< C—1| VwoIIl 2 ©! • exp{C0 2 }. 

[I 


M 4 p 8 M 2 M|p 10 ^ p 4 Mfp 4 „ r2 , 2 

—T- + + P, ©2 = ~r T —+ M 2 M 4 p 2 . 

p 6 p 6 — p* p- 


Proof. Take t-derivative of the momentum equation, 

(3.31) pu tt + (pu)-Vu t -div (2p(p)d t ) + S7 P t = -p t u t -(pu) t -Vu + div (2 p(p) t d). 

Multiplying (I3.31j) by tu t and integrating over 12, we get after integration by parts 
that 


(3.32) 


— J p\u t \ 2 dx + 2t J p(p)\d t \ 2 dx 
= —t Pi\ u t\ 2 dx — t ( pu) t ■ Vw ■ u t dx — t 


2p(p)td ■ Vu t dx . 
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Let us estimate the terms on the righthand of (13.321) . First, utilizing the mass 
equation and Poincare’s inequality, one has 


(3.33) 


-t J p t \u t \ 2 dx 
= —2 1 J pu ■ V u t ■ u t dx 

< Op • t ■ \\^pu t \\ L 3 ■ \\Vu t \\ L 2 ■ \\u\\ L 6 

< C-p * • t ■ \\y/pu t \\l 2 ■ \\y/pu t \\l 6 ■ ||VUt || L 2 • ||Vu || L 2 

< C~p^ ■ t ■ \\y/pu t \\ 2 L2 • ||Viitllla • ||Vn|| L 2 

< ^p-t\\Wu t \\ 2 L 2 +Cp~ 3 -t\\^pu t \\l 2 - ||Vu||la. 


Second, utilizing the renormalized mass equation for p(p), 
— 2 tj p(p)t ■ d ■ S7u t dx 

(3.34) < °t J \u\ • |V/i(p)| • |d| • |Vu t | dx 

< ct ■ ||V/i(p)|| L 3 ■ ||Vni || L 2 • \\d\\ L 6 • ||n|| L o 

< CMt ■ \\Vu t \\ L 2 ■ \\Vu\\ 2 m . 

It follows from Lemma 13.41 that 


(3.35) 


— 2 1J p(p)t ■ d ■ \/u t dx 
1 C! M 2 

< -p ■ t\\S/u t f L 2 + — ■ t {M*p 2 \\^pu t \\ 4 L 2 + Mlp 8 \\Vu\\f 2 ) 


p, 


< -p-t\\Wu t \\ 2 L 2 + 


CM 2 M 4 p 2 

p 


■tWy/PkWh + 


CM 2 Mlp s 

p 


■ t|| Vu 


112 
I L 2 - 


Finally, taking into account the mass equation again, we arrive at 
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Hence, it follows from Sobolev embedding inequality, Gagliardo-Nirenberg in¬ 
equality, and Lemma 13.41 that 

(3.37) 

J\ < Cp-t\\u t \\ L 6 ■ \\u\\ L 6 • ||Vtt||| 3 
<Cp-t\\Vu t \\ L 2 -\\Vu\\l 2 -\\Vu\\ H i 

< Cp-t\\Vu t \\ L 2 ■ ||Vu|| 2 2 ■ (M 2 p?\\y/pu t \\ L 2 + M 2 p 2 ||Vu||| 2 ) 

< L • i||v« t ||| a + GUd. i|iv^tlli= • ||v«iii, + —• i||v«nll 

o~ fl 

Similarly, it holds that 

(3.38) 

J 2 < Cp ■ t\\ut\\ l 6 ■ || x, 2 ■ \\u\\ 2 L e 

< Cp-t\\Vu t \\ L 2 • ||Vtt|| H i • ||V-u || 2 2 

1 Il9 CMfp 3 .. ^ ll9 

< -fl ■ t||VUt|| L a + -— • tWy/pUtWtf 

o~ fl 

and 

(3.39) 

J 3 < Cp-t\\Vut\\i^ ■ l|Vw|| L 6 • ||m|| 2 6 

< L-iiiv«,iii, + CSeL . t \\^pu t \\i, . iiv«h | 2 + £*%£ ■ *iv«ig 

O f-L jh 

Owing to Lemma 13.41 and Sobolev embedding inequality, 

J± < Ct\\y/pu t \\ 2 L 4 • ||Vn || L 2 

(3 40 ) < CtWy/putWli-ptWVutWl, ■ \\Vu \\ L 2 

< L ■ «irvein. + ^«■vuiiu 

o — fl 


•||Vn||i 2 + 


4-6 


CM$p { 


■ t|| Vm 


110 
I L 2 > 


Combine all the above estimates fl3.33MI3.40lb 


(3.41) 


J p\u t \ 2 dx + p- t\\Vu t \\ 2 L 2 

< • t\\y/ptH\\b ■ IIV^HI, + - • t\\^u t \\ 2 L 2 • ||Vu|| 4 i2 


p 3 


CM£p 6 


p 


tL 

,10 , CM 2 M$p 2 
IL 2 


H--— ' i||VM||i u 2 H • t||\/p M t|| 4 

p 

CM' 2 M^jf 


L 2 


+ 


p 


■ *||Vu||# + / p\u t \ dx. 
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Applying Gronwall’s inequality, 


sup t / p\u t \ 2 dx + p / t\\S7u t \\\ 2 dt 

te[o,T] J —Jo 

T / CM 2 p fs 


< 


L</o \ P 

t r 

• exp 


■t ||Vu||i u 2 + 


10 , +IIT7 .. i|12 , „ l|2 


P 


•t||V'u|| L 2 + Wy/puth* dt 


Cp* CM$p 3 


/A 


+ 


p 


4 CM 2 M A p 2 


l|V«||i 2 + 


P 


IIVMII 


L 2 


dt 


Taking (I3.2(jp and (I3.23P into account, 


sup t / p\u t \ 2 dx + p / t\\Wu t \\ 2 L 2 

te[o,T] J Jo 


(3.42) < 


T (C-Mffi Ilin CM 2 Mlp s Ill2 _ ll2 . , 
*||Vu|| i2 H y-• ^||Vm|| L 2 + \\yfpu t \\ L 2 ) dt 


K 




p a p 2 


According to Theorems 13.211X51 and the assumption 113 .1 8H . 


,T pT 

t\\Vu\\)®2 dt < sup t||V?i ||^2 • sup \\Vu \\ 6 L 2 • / \\S7u\\ 2 L 2dt 


< 


te[0,T] 

C-p 2 


te[o,T] 


- 7;2 


P 2 


Kill* < ^-AllVfolH.. 




Similarly, 


*||Vix||^<^^||V« 0 ||l 2 

P 


And by virture of Theorem 13.31 


\\VpUt\\ 2 L 2 dt < p ■ \\Vu 0 \\ 2 l2 . 


Hence, 

(3.43) 


sup t / p\u t \ 2 dx + p / t\\S7u t \\ 2 L 2 dt 

te[o,T] J —Jo 


p- 


p- 


^ I|VT ll2 /M 2 4 p» M 2 M»p w \ f fjf M 2 p A 2 4 _ 2 

< G||Vmo||l 2 ( —yy—I ~y-h /i ) • exp C ( — H —-h M M 2 p 


p A p 2 
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On the other hand, multiplying (13.41)1 by t, one has 


(3.44) 


j t Cj f pWt\ 2 dx) + pt 2 \\Wu t \\ 2 L 2 


Cp 3 


2t;3 


n m 2 7> 

< ^ ■ GlVPX.ll!> ■ l|Vx||!> + —■ GlVP“<ll!> ' HVulfc 


p 3 


p 


+ 


+ 


CMf-fi 


*2 || V7„. II10 , 1 2II r^„. 114 


P 


■ t ||VnUyi + 


p 


■nvputWh 


CM 2 M$p 8 +2||V7 ||12 , , , |2 


p 


■t ||Vu ||^2 +t p\u t \ dx. 


Applying Gronwall’s inequality, 

( 3 - 45 ) 

sup t 2 / p\u t \ 2 dx + p / t 2 \\S7u t \\ 2 L 2 dt 
te[o,T] J —Jo 

T (cm?? . t 2 nv«ni» 3 + OAMld. i 2 nv«iiB + *n vp«<III>) dt 


< 


IJ o 


p 


■ {C [ d + MM. + M‘Mt? 

.. P P 


P 

2 n ,ri-X2 


According to Theorems 13.21 and 13.51 


t 2 || Vw ||}°2 dt < sup t 2 ||V«||l 2 • sup ||Vn || 4 L 2 • / ||V-u|| 2 2 dt 


< 


te[o,T] 

Cp 3 


te[o,T] 


P d 


|| Vu 


2 

Oil L 2 ' 


Similarly, 


Hence, 

(3.46) 


t 2 ||V M ||i 2 2dt<^||V Mo ||i 2 

P 3 


sup t 2 / p\u t \ 2 dx + p / t 2 \\S7u t \\ 2 L 2 dt 

i 2 g [ o , t ] J —Jo 

r2 AjS-Ttll 


< C\\Vu 0 \\ 2 l2 


2 , , m 2 m$p 11 , \ _ , m 2 V , 2 


p 4 


+ 


p 4 


+ T exp fU + 7 


+ M z Mfp z 


□ 


Lemma 3.7. Suppose (p,u,P) is the unique local strong solution to (11.11) on 
[0,T], with the initial data (p 0 ,no), and satisfies the assumptions (13. 17jl - (13. 1911 . 
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Then for any r e (3, min{g, 6}) 
(3.47) 

[ ||V«|| L °° dt 


< C||Vmo||l2 


_ \ I 6(r-l) 

5r_6 3(r—2) / Q \ 2 I „ ^ 5r-6 Q r 

M r p 4r /i ^ ( 1 H— ) 0? exp{(709} + M r r - 

“ /V T 


Proof. By virture of Lemma [2.11 one has for r G (3, min{g, 6}) 


HVuHiyi.r < CM r (\\pu t \\ L r + ||/ra • V-uHir) 


(3.48) 


3 (y 2) 


< CM r \\pu t \\ L 2 r • llpMtll L 6 2r + p\\u\\lfi • || Vn 11 L 6r/(6-r 


3 (y 2) 


r-1) 


<cm, iip« i |i 1 ? +p|iv«n 1 ! r« • uvular,: 


Applying Young’s inequality and Sobolev embedding inequality, 

(3.49) 

Sr — 6 6 — r 3(r-2) 5r-6 5r ,_ 6 6(r-l) 

\\X7u\\ w i,r < CM r p ir ■ Wy/putWtf • ||Vut|| i2 2r + CM r r p T • ||Vn|| i2 r , 


Hence, 


|| Vu||l°° dt 


(3.50) <C ||Vu|| w i,rdi 


Jo 

f T ( 5r —6 SziT 3(r —2) 5r-6 5r _ 6 6(r-l) ’ 

< C J [M r p 4r \\y/put\\l 2 II V'U t || L 2 2r + Mr P r 1111 ^ 2 " | dt. 


If T < 1, according to Theorem 13.61 


3(r —2) 


WVpMv • l|VM t || i 2 2r dt 


< 


* 2 II\/MIU 2 


b —r 
2 r 


3(r-2) 

C || V'Ut||i2 ) “ r ■ t~3 dt 


6—r 
2 r 


3(y—2) 


< C [ sup ^Hv/p^ll^ 
ye[o,T] 


*l|Vut||| 2 dt 


< Cp-^ 1 • HVitollia©? exp{C0 2 }. 


rT 2 r N ^ 

t h-6 dt 
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If T > 1, applying Theorem 13.61 again. 


3Q--2) 


\WpM£ •||V^|| i 2 2 ’' dt 


< 


o—r 
2 r 


t 2 \\VpU t \\ L 2) • f2||V« t || L 2 -t *dt 


3(r-2) 

2r 


3(r-2) 


+ / (t||\/MllL 2 ) 2r ■ (t||Vu t || L 2) ■t 1 dt 


6 —r 
2 r 


< C [ Slip t*\\y/pUt\\L* 
ye[o,T] 


C ( sup t\\y/pu t \\ L 2 
te[o,T] 


t\\Vu t \\ 2 L 2 dt 


3(r—2) 
4 r 


o—r 
2 r 


t 2 || dt 


< Cfi-^WVuoWvGf exp{C© 2 } + CpT^ ( P - ) ||Viio|U 2 ©? exp{C0 2 } 


d 


3(r-2) 

4r 


t r + e dt 


m 6 + 7- 

r T \ 4r- 

t r + 6 dt 


'i 


< Cp 3( 4 r } ( 1 + ^ ) I|V«o||l 2 © 2 exp{C'© 2 }. 


On the other hand, 


(3.51) 

Therefore, 

(3.52) 


6(t—1) 


||Vn|| i2 r dt < sup \\Vu\\j 

te[o,T] 


||V«||| 2( it<^||V«o||L2 

/i 


/ || Vm||l°° dt 
Jo 

< C'H Vt/o||z,2 

±C 2 (M,p,p,p)\\Vu 0 \\ L 2 . 


5r—6 3(i—2) 


_ x 1 1 6Q—1) 

M r p p Hr- I 1 + — ) Of expiC*©?} + M r r — - 

d) d 


Theorem 3.8. Suppose (p,u,P) is the unique local strong solution to (II. Hi on 
[0,T] and the assumptions (j3.17p - fl3.19]) . There exists a positive number eo, de¬ 
pending on Q, q, M, p, p, p such that if 



1 Vrto| L 2 < O , 

then 


(3.53) 

sup ||V/i(p)||l9 < 2 M , 
te[o,T] 

and 

sup || Vp||z,« < 2|| V/9 0 |U^ 

ie[0,T] 
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Proof. Consider the Xj-derivative of the renormalized mass equation for p(p), 
(difj,(p)) t + (diU ■ V)p(p) + u ■ Vdip(p) = 0 . 

It implies that for every t G [0, T\, in view of Lemma [3.71 one has 

^ 3 54 ) l|V/z(p)(t)|U« < ||V/i(po)||m -expj^ ||VM|| L oods| 

< ||V//(p 0 )IU« ' exp {C 2 (M, p, p,p) • || Vitollr 2 } ■ 
Choose some small positive constant eo, which satisfies 

e 0 < i ejj ■ CipT 2 (M| + M 4 M 2 6 )p 4 < In 2, 

and 

e 0 • C 2 (M,p,p,p,) < In 2. 

If ||Vn 0 || L 2 < e 0 , (13.53P holds. 

Similarly, 


(3.55) 


||Vp(£)||l? < ||Vp 0 |U« - exp | f || V«||l°° ds 


< 2||Vpo||L9, 

Therefore. Theorem 13.81 is proved. 


3.2. Proof of Theorem 11.21 With the a priori estimates in Subsection 3.1 in 
hand, we are prepared for the proof of Theorem 11.21 

Proof. According to Theorem 11.11 there exists a T, > 0 such that the density- 
dependent Navier-Stokes system (II.ip has a unique local strong solution (p, u, P ) 
on [0, T*]. We plan to extend the local solution to a global one. 

Since ||V/x(po)||z,« = M < AM, and due to the continuity of Vp(p) in L q and 
Vuo in L 2 , there exists a T\ G (0,T*) such that sup 0</<T] ||V//(p)(t)||L« < AM, 
and at the same time sup 0<t<Tl ||Vm(£)||z ,2 < 21| Vrzo||z. 2 - 
Set 

T* = sup {T| (p, u, P) is a strong solution to (11.11) on [0, T]} . 

(p, u, P ) is a strong solution to (11.ip on [0, T], 
sup \\X7p(p)\\ L q < AM, and sup ||Vu||x ,2 < 2||Vu 0 ||l 2 - 

o <t<T o <t<T 

Then T* > 7\ > 0. Recalling Theorems 13.31 and 13.81 it’s easy to verify 

(3.56) T* = T*, 

provided that ||Vmo||l 2 < eo as assumed. 

We claim that T* = oo. Otherwise, assuming that T* < oo. By virture of 
Theorems 13.31 and 13.81 for every t G [0, T*), it holds that 

(3.57) ||Vp(f)||i? < 2||Vp 0 ||i8, and ||Vu(f )|| L 2 < a/ 2|| Vw 0 || L 2 . 

which contradicts to the blowup criterion (j 1.5 p . Hence we complete the proof for 
Theorem 11.21 



T* = sup < T 


□ 
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